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accepted for publication in the Journal. 

11 We quickly print-on-demand an issue every 3 papers, in 
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12 We do not demand from the author to transfer his own copy- 
right to us. 

13 Nevertheless, we publish only an original research work 
paper and only if the author does assure us that he owns 
the copyright of his own paper and submit this paper or a 
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even for republication under the Creative Commons-Attribu- 
tion 3.0 Unported License, CC-BY 3. 
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The Very First Ever Made 
3D/4D Virtual Laboratory for Antennas 

N.I. Yannopoulou, P.E. Zimourtopoulos * 

Antennas Research Group, Austria - Hellas [l, 2] 

EECE Dept, Democritus University of Thrace, Hellas [2] 

Abstract 

Based on the experience we have gained so far, as indepen- 
dent reviewers, we thought that may be proved useful to pub- 
licly share with the interested author, especially the young 
one, some practical implementations of our ideas for the in- 
teractive representation of data using 3D/4D movement and ani- 
mation, in an attempt to motivate and support him in the de- 
velopment of similar dynamic presentations, when he is look- 
ing for a way to locate the stronger aspects of his research 
results in order to prepare a clear, most appropriate for 
publication, static presentation figure. For this purpose, we 
selected to demonstrate a number of presentations, from the 
simplest to the most complicated, concerning well-known an- 
tenna issues with rather hard to imagine details, as it hap- 
pens perhaps in cases involving Spherical Coordinates and Po- 
larization, which we created to enrich the very first ever 
made Virtual Laboratories of Antennas, that we distribute 
over the Open Internet through our website Virtual Antennas. 
These presentations were developed in a general way, without 
using antenna simulators, to handle output text and image 
data from third-party CAS Computer Algebra Systems, such as 
the Mathematica commercial software we use or the Maxima 
FLOSS we track its evolution. 



Introduction 

The very first ever made 
Virtual Laboratories of An- 
tennas, which are available 
on the Internet through Vir- 
tual Antennas website [1], 
are exclusively based on our 
alternative form of Antenna 
Theory [2], and are founded, 
by applying the learning by 



teaching method to Antennas 
education [3], in order to 
support the quick comprehen- 
sion of 3D space matters un- 
avoidably related to anten- 
nas. The initial material of 
Virtual Antennas was prepared 
during the years 1996-1997 
and presented to Antenna stu- 
dents during the spring seme- 



ster of 1998 [4] . On 1999, 

the existence of this visual 
material was announced to the 
first ever appeared on the 
web EMLIB Electromagnetics 
Library, maintained those 
years by Jet Propulsion Labo- 
ratory of NASA [5] . On the 
same year, all of the mate- 
rial was approved for inclu- 
sion to MathSource repository 
of Mathematica [6] . 

Since 2000 and until to- 
day, the Virtual Antennas ma- 
terial is constantly in- 
creased and improved. More- 
over, a number of websites 
either use this material or 
suggest its use [7] . On 2009, 
a sample of its current de- 
velopment state was approved 
for publication by the Wol- 
fram Demonstrations Project 
[8], an event that was an- 
nounced in our Creative Com- 
mons Network web pages [9] . 

During the last years, our 
voluntary reviewing work re- 
vealed some needs of poten- 
tial authors -especially the 
young ones- related to the 
presentation of their re- 
search results, which perhaps 
may be confronted by similar 
to our interactive presenta- 
tion techniques, and thus we 
decided to present in this 
paper our continuing work for 
the Virtual Laboratories of 
Antennas . 



Presentation Development 

and Use 

The implementation of a 
presentation idea demands the 
expression of the theoretical 
idea formulation in the CAS 
language we select to produce 
output data in the appropri- 
ate file format for movement 
or/and animation (we would 
like to emphasize that in the 
case of Mathematica we use, 
while we abandoned the pack- 
age written by Novak [10], we 
still run the one written by 
Donley [11] ) . These files 
are: (a) WRL (WoRLd: Open by 

Web3D Consortium), a plain 
text file with a known struc- 
ture described in the VRML 
Virtual Reality Modeling Lan- 
guage [12] for a visual ob- 
ject, (b) a bunch of images 
in a selected non-destructive 
format; usually BMP, and (c) 
NBP text file (Note Book 
Player: Mathematica propri- 

etary). After that, the de- 
velopment of the presentation 
requires: (a) either a text 

editor to correct bugs or to 
add special features into WRL 
output file, (b) a graphics 
editor to align images or 
correct blemishes, as well 
as, a video editor to handle 
these images as frames of the 
final movie file in a se- 
lected format, usually AVI 
(Audio Visual Interleave: 
Microsoft proprietary) . It is 
worth to notice that accord- 
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ing to our experience, the 
non-computing time needed: 
(a) to design the idea imple- 
mentation, has been reduced 
by the years, from a couple 
of weeks at the beginning to 
just a couple of days now, 
and (b) to develop the pre- 
sentation, after a tedious, 
routine work of editing each 
frame, varies considerably 
according to the frame theme, 
e.g. from 1.5 minutes for 
that in Fig. 7, to 10 minutes 
of Fig. 10. 

The software needed for 
presentations to work under 
MS Windows is: a web browser 
(MS -IE 3 . 02r+) , a VRML viewer 
add-on to that browser (World- 
View 2.1), an AVI player 
(Mplayer2), and the Mathema- 
tica 7 Notebook Player. 

3D/4D Presentation Samples 

In order to abbreviate the 
figure descriptions in the 
following, we have to notice 
that: (a) for every animation 
presentation, a typical full- 
window screen capture of the 
recommended AVI player is 
shown after its pause button 
was pressed, (b) for every 
movement presentation, a ty- 
pical full-window screen ca- 
pture of the recommended on- 
line WRL viewer, with a cur- 
sor perhaps to indicate the 
existence of an additional 
feature, and (c) for both 
presentation types, the basic 
colors (R, G, B) are used in 



that order to correspond ei- 
ther to the CCS Cartesian Co- 
ordinate System (x, y, z) 
axes, (xoy, yoz, zox) planes, 
and its unit vectors, or to 
the SCS Spherical Coordinate 
System (r-radius, theta-semi- 
circle, phi-circle) curves, 
(r-sphere, theta-semi-cone, 
phi-semi-plane) surfaces, and 
its unit vectors, while the 
same color correspondence 
holds for the radiation pat- 
tern cuts by the mentioned 
CCS planes or SCS surfaces. 

After that, short descrip- 
tions of the samples, are 
following, while all of these 
presentations will be always 
available in authors' group 
repositories in Virtual An- 
tennas [13] and Google Code 
Repository [14] . 

Fig. 1 shows 1 frame out 
of 12 of an AVI repeated-for- 
ever animation for the 3D 
plane-time presentation of 
the considered as time-har- 
monic sinusoidal current wave 
amplitude, on a 3 wavelengths 
portion of a long thin wire 
loop antenna terminated on 
some complex impedance. The 
propagated current wave p is 
decomposed in the dual cou- 
ples: (incident wave i, re- 
flected wave r) and (standing 
wave s, transmitted wave t), 
while the imposed letters and 
arrows on the figure indicate 
these current waves as well 
as the direction of their mo- 
tion . 



N.l. 

Fig. 2 shows a screen-shot 
of a WRL presentation with 5 
additional predefined view 
points of 12 frames anima- 
tion, as the cursor indica- 
tes, of a 4D space-time asyn- 
chronous movement, from a 
view-point in the lst-octant. 

In essence, this presentation 
is an alternative geometric 
meaning of the considered as 
standing wave current on an 
adequate length of a thin 
wire dipole antenna. 

Fig. 3 shows a screen-shot 
of a WRL presentation of 4D 
space-time asynchronous move- 
ment from a view point in the 
lst-octant, of the 3D differ- 
ence volume element formed by 
adjacent SCS coordinate sur- 
faces and lying between the 
1st and 2nd octant. 

Fig. 4 shows a screen-shot 
of a WRL presentation of a 4D 
space-time asynchronous move- 
ment from a view-point in the 
lst-octant, for the triples 
of SCS unit vectors in the 
shown directions, that is 13 
triples on the CCS main- 
planes, as well as, 2 unde- 
fined triples on the irregu- 
lar z-axis of the CCS-to-SCS 
transformation . 

Fig. 5 shows a screen-shot 
of a WRL presentation of a 4D 
space-time asynchronous move- 
ment from a view-point in the 
lst-octant, for all 3 SCS 
surfaces, curves, and unit 
vectors in a direction of the 
1st- octant. 
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Fig. 6 shows 1 frame out 
of 37 of an AVI animation for 
the 4D plane-time presenta- 
tion of a SCS phi-circle cur- 
ve as the intersection bet- 
ween 2 SCS surfaces, that is 
of a r-sphere and a theta- 
semi-cone, while theta angle 
-which changes from 0° to 
180°- has the value of 50°. 
Different colors illustrate 
the inner and outer surfaces 
of the r-sphere and a cut off 
spherical section between ad- 
jacent meridians permit us to 
see the interior of that 
sphere . 

Fig. 7 shows a screen-shot 
of a WRL presentation of a 4D 
space-time asynchronous move- 
ment for three different po- 
larizations at three differ- 
ent points: linear, circular 
and elliptical on each of the 
coordinate axes x, y and z 
respectively to emphasize 
that, in general, the pola- 
rization of an antenna is not 
constant . 

Fig. 8 shows 1 frame out 
of 20 of an AVI animation for 
the 4D space-time presenta- 
tion of the non-linearly CCW 
polarized time-harmonic real 
electric field E, from an an- 
tenna, as well as, its decom- 
position to two linearly po- 
larized time-harmonic real 
electromagnetic fields, as 
they are uniquely defined by 
the constant of time linearly 
independent space vectors E c 
and E s . 
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Fig. 1: Current waves 



Fig. 4: SCS: Unit vectors 









Fig. 2: Standing current- 
voltage waves relation 
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Fig. 3: SCS: Volume element 



Fig. 5: Coordinate curves, 
surfaces and unit vectors 




Fig. 6: Coordinate surfaces 
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Fig. 9 shows the last fra- 
me out of 23 of an AVI anima- 
tion for the 4D space-time 
presentation of the path 
traced by the electric far- 
field arrow tip, which is CW 
elliptically polarized in the 
direction of y-axis. 




Fig. 7: Polarization: 
Directional dependence 
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Fig. 8: Polarization: 
Non-linear CCW 



Fig. 10 shows 1 frame out 
of 19 of an AVI for the 4D 
space-time presentation of a 
normalized CCW elliptically 
polarized electromagnetic far 
-field in the shown direction 
of propagation, from a point 
of view in the 1st octant of 
an observer mirrored to (r, 
theta) plane. 




Fig. 9: Polarization: 
CW elliptical 




Fig. 10: Polarization: 
Far-Field CCW elliptical. 
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Fig. 11 shows 3 frames out 
of 73 of an AVI for the 4D 
space-time presentation of 
the direction dependency of 
antenna polarization from a 
center fed short crossed 
dipoles on y and z axes, 
which have a current phase 
difference of 90° . The three 
different polarizations are 
shown with their orientation 
for an observer on the xOy 
plane (theta = 90° and in 

three different phi angles of 
0°, 240°and 270° . 

Fig. 12 shows 4 frames out 
38 of an AVI animation for 
the 4D space-time presenta- 
tion of the space radiation 
pattern and its three main- 
planes cuts, for the specific 
case of an antenna consisting 
of 2 dipoles, with the same 
direction of the unit vector 
(0.2, 0.4 0.894), each of 

2.4A length long, with their 

centers placed 0.25A apart on 
an axis with a unit direc- 
tional vector (0.3 0.5 

0.812), fed with a 30° cur- 
rent phase difference. 

Fig. 13 shows a screen- 
shot of a WRL presentation 
with 1 additional predefined 
view point of 12 frames ani- 
mation possibility, as the 
cursor indicates, of a 4D 

space-time asynchronous move- 
ment, from a view-point in 
the lst-octant. This presen- 
tation concerns the interior 
of an anechoic chamber and 

imitates a discone antenna 



rotation in order to visual- 
ize the way by which the mea- 
surement of a plane-cut of 
its radiation pattern is ac- 
complished . 
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Fig. 11: Polarization: CW 
Circular - CCW Elliptical - 

Linear 
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Fig. 14 shows a sophisti- 
cated combination of a move- 
ment with an animation pre- 
sentation, which in fact is 
an application that was de- 
veloped using the version 7 
of Mathematica, in a way that 
permits the definition of any 
(theta, phi) direction in 
space for a dipole of vari- 
able length. Notably, this 
application has all the avai- 
lable VRML movement features, 
with 6 predefined view 
points, while 3 different an- 
imations can be run simulta- 
neously. In addition, the 



mesh, opacity, rendering goal, 
and the evaluation step for 
the radiation pattern can be 
also defined. 

Fig. 15 shows the last 
frame out of 100 of an AVI 
for the 4D space-time presen- 
tation of a detailed study 
for some of the linear dipole 
characteristics regarding the 
plane and space radiation 
patterns, the direction from 
the dipole axis of their 
first common maximum, the di- 
rectivity, and the input ra- 
diation resistance; this re- 
ally is a 4-tuple of frames. 




Fig. 12: Radiation Pattern: 3D and its 2D main-plane cuts 
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Fig. 13: Anechoic chamber - discone antenna 
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Fig. 14 Radiation Pattern: A composite 3D/4D presentation 
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Centre -Fed Linear Dipole with Standing Wave of Current 



Normalized E Radiation Pattern 




ft . Plane Pattern B . Space Pattern 
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Fig. 15: Antenna characteristics 



Conclusion 

The Antennas Research Group 
(ARG) carries on for the en- 
hancement of this low-cost 
Virtual Laboratory with new 
material placed to web site 
www.antennas.gr [13], as well 
as to Google Code Repository 
[14], and Fig. 13 and Fig. 14 
are only two of many more 
possible examples. 

The evolution on computer 
systems in recent years has 
increased the demand for new 
teaching, learning and re- 



searching methodologies, which 
this paper supports by de- 
scribing available material 
for a web-based Virtual Labo- 
ratory for Antennas that pro- 
vides interactive computer 
aided visualization in vari- 
ous Antenna Theory topics. 

The main goal is con- 
stantly attained since this 
continuously developed mate- 
rial cannot be produced by 
the means of the conventional 
literature . 
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Abstract 



We have developed two methods to study an inverse 
scattering problem of electromagnetic waves in chiral media, 
for a spherical scatterer perfect conductor. The first 
called Far-field inverse problem was described analytical in 
our previous work. The second called Near-field method is 
given in detail in the present paper. It is a geometrical 
method based on the scattered fields. Using Near-field 
experiments, in which the scattered field is measured at the 
source, we solve the corresponding inverse scattering 
problem that is to determine the coordinates of the center 



and the radius of the spherical 

Introduction 

We study an inverse scat- 
tering problem of electromag- 
netic waves with time harmo- 
nic dependence applying the 
method of Near-field [l]-[3]. 
Spherical electromagnetic wa- 
ves generated by a point 
source incident on a spheri- 
cal scatterer perfect condu- 
ctor in a chiral media. 

Knowing the incident and 
scattered wave fields in the 
inverse problem we are 
looking for the coordinates 
of the sphere center and its 
radius . 

The corresponding problem 
is solved by geometrical me- 
thods using either energy 
scattering cross-section, or 
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scatterer . 

the scattered fields. In the 
first case we refer to the 
inverse Far-field problem [4] 
and in the second case to a 
Near-field inverse problem. 

The chiral materials exhi- 
bit the phenomenon of optical 
activity i.e., the phenomenon 
that in which the plane of 
polarizations of linearly po- 
larized light is rotated as 
the light passes through an 
optically active medium. 

The chirality is a proper- 
ty which is often found in 
nature and reflects the asym- 
metry in the spatial inver- 
sion. An electromagnetic wave 
into a chiral media is analy- 
zed in a counter-clockwise 
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( LCP ) and a clockwise (RCP) 
Beltrami field. Using vector 
spherical harmonic functions 
and appropriate expansions of 
the Beltrami fields and by 
extension the development of 
spherical electromagnetic wa- 
ves in spherical wave fun- 
ctions, we calculate the 
exact solution of the scat- 
tering problem for a spheri- 
cal perfect conductor in chi- 
ral media, and the correspon- 
ding back scattering [3] -[4]. 

In the inverse problem we 
calculate the scattered Near 
-field, when the source is in 

position r = r 0 . A similar pro- 
blem for achiral materials 
has been studied in the works 
[ 1 ], [ 2 ]. 

Statement of the problem 

Consider a point source at 
the position r 0 that produces 

spherical electromagnetic wa- 
ves in a chiral media near a 

scatterer Q ", perfect con- 
ductor, i.e. the surface of 
the boundary condition is sa- 
tisfied 



A spherical incident ele- 
ctromagnetic wave (E“ c (r), 

H“ c (r)) with time harmonic de- 
pendence in accordance with 
the Bohren transformation, 
analyzed in spherical Beltra- 
mi fields Q^(r) and (r) 
as follows: 



E JDC (r) = 

fo 


Q^(r) + Q“r 0 ( r ) 






< 


1 






(2) 


C(r) = 


— (Q , n r O') - 


-Qr™ 


;(r)) 


V 


lq 0 


r i 






where q 


= (P / £ ) 1/2 


is 


the 


in- 


trinsic 


impedance 


of 


the 


chi- 



ral medium. The Beltrami 
fields satisfy the equations, 
[ 5 ], [ 6 ], 

VxQ 4re (r) = y,Q 4r0 (r) 

VxQ Rr0 (r) = -y R Q 4rs (r) 



where y L , y R are wave numbers 

for Beltrami fields and are 
given by, 



k k 

Yl ~ 1 - kp ' Yr ~ 1 + Arp 




n x (E i,,c (r) + E sc (r)) = 0 

( 1 ) 

r = a 

We consider that the exte- 
rior Q + = R 3 /Q" of the scat- 
terer is homogeneous chiral 
media with fixed chirality 

(3, dielectric constant e and 
magnetic permeability p . 



With k = co(e p) 1/2 , co being the 

angular frequency. The indi- 
ces L and R denote the LCP 
and RCP fields respectively. 
The spherical incident Bel- 
trami fields with suitable 
normalization have the form, 
as defined in the following 
works: [4] (issue 2, p. 9, 

relations (4), (5)), and [3]. 
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If (r) and Q“‘(r) are 

Beltrami fields corresponding 
via transformation to Bohren 

in Ef c and Ht nc , then the 

scattering problem for the 
perfect conductor can be 
formulated in the following 

way: be found Q'“'(r), Q“‘(r) , 

which belong to the space 

C 1 (Q + ) n C(Q + ) , such that: 



(E* c (r), H* e c (r)) . The scattered 
electric field will be depen- 
ded on the polarizations , 

p R , (see [4] (issue 2, p. 

9, relation (6)) and will have 
the decomposition 

(r | p LI p R ) = Q“ r (r | p L , p R ) + 

( 6 ) 

+ Q£r. ( r I Pi./ P R ) 



(i): 



v x Q™ (r) = y L q[% (r) 

V x Q™ (r) = -y R Q£.(r) 



r e Q 



(ii) : n x Q^(r) =-n x Q^(r) 

r e S =dQ~ 



(5) 



f x Q£. (r) + iQf ( c r< (r) = o(J) 
r x Q» c r „ (r) - iQf, (r) = o(h 



r -> oo 



The limits on radiation 
conditions (5(iii)), are ta- 
ken uniformly in all dire- 
ctions reS 2 , where S 2 is 

the unit sphere in R 3 and n 
is the outward normal unit 
vector perpendicular to the 
surface on the scatterer. 

The incident electromagne- 
tic wave (Ef" c (r) , H“ c (r)) on the 

scatterer Q~ generates the 
corresponding scattered field 



Where Q“ e (r|p„p R ) and 

Q£ r . ( r I Pu P r ) are the corre- 
sponding scattered Beltrami 
fields which have the fol- 
lowing behavior, when r -» oo , 

[5], [7], 

Q" ro (r I P„ P R ) = K (Y A r)- 

i {7) 

■9A, r Ar\p L ,p R ) + 0(±r) 

with A = L, R . The functions 

g, . and g„ r are the LCP and 

RCP far-field patterns respe- 
ctively, which are defined by 
the following relationship 
[7], 

g^,r e (r|p„P R ) = 

{fix [y A V X E® c (r' I P L , p R ) - (8) 

S 

- m A y 2 E^(r' | p„ p R )e-^ r r c/s(r') 

-1, A = L 

with m A ={ , y 2 = y L y R . 

1, A = R 

V 
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Near-field Inverse problem 

Using the exact solution 
found in the paper [4] (issue 
2, pp. 10-12) and considering 
appropriate asymptotic forms 
of the Bessel and Hankel 
functions in low frequencies, 

i.e. y A a «: 1 , calculate asym- 
ptotic expressions of the 
coefficients a L n , a* , b L n and 

b R „ . Specifically we have [3], 



i (1 + (3/c) n + 2 2n + , 

2 nC n (l-p/O n + 1 ’ 



l l + P k , . 

a„ o — (y,a) 

n 2ifi(2n + iy rL 

y.a -> 0 



2n + 1 



(9) 



i (l-p/cr 2 2n + 1 
2 nC„ (1 + P k) n + lKYi ’ 



bt ~ - 



y R a — > 0 



( 10 ) 



i (1 - PfO 
202 n + 1) 



(/ R a) 



2n + 1 



Where 

£„ = l-3-5---(2 n - 1) = (2n)!/(2"n!) 

In the inverse problem we 
calculate the scattered Near- 
field in the source r = r 0 . By 
the relation (21), in work 
[4], with r = r 0 and for LCP 
incidence, we obtain [3], 



E* e c (r 0 |p L ) = 

= J n(n 2 + 1) B L n a L n {(h n (y L r 0 )x + h n (y L r d ) qj) + i(h n (y L r 0 )x - h n (y L r e ) iji) } + 



( 11 ) 



+ J n(n 2 + ^ B n a ni( h n(Y R r B )x ~ h n (y R r e ) iji) + i{-h n (y R r 0 )x - h n {y R r d ) iji)} 

n= 1 



Using the asymptotic rela 
tions : 

MVal 

Mv A r B ) ( y A rJ 



IUYa^} zHtn 
h B (Y A r B ) ( y A r e ) n + 1 

h n (y A a) ~ - y 2 

i ( y A a ) 



with y A r B -» 0 , A = L, R and the 
series, 



.2 n + 1 



T 

1-T 2 



( 12 ) 



2/1 + 1 



r 3 (l + t 2 ) 
(1 - t 2 ) 3 



(13) 



+ (14) 
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Y n(2n + l)r 



2/1 + 1 



r(3r 2 + r 4 ) 
(1 - r 2 ) 3 



^ (15) 



2n + 1 



Z £-l I 1 

n 



.2 n + 1 



(16) 



(1 - TO^ft/c) T 5 

2(y A a) 2 (1 - t 2 )' 

y aO — ^ 0^ >A — i—f r 



( 20 ) 



_ r[r 2 In (1 - t 2 ) + 2r 2 - In (1 - r 2 )] 

“ r 2 - 1 

and because of the relation- 
ships : 



1 - Pft 
i + 3^ 



, 0 < r < 1 , 



where r = a/r 0 , we obtain from 

( 11 ) [ 3 ], 



i pj ~ 

(i + gft) r 5 

2(/ t a) 2 (1 - t 2 ) 3 ' 



(17) 



„ (1 + B/c) r 5 

r ' UIPJ 2 (y L a) 2 (1 - r 2 ) 3 

y L a -> 0 



(18) 



Similarly for the RCP inci- 
dence, we have that 



E sc (r I p ) ~ — — ^ — L 

r 0 \ 0 I P/? / 2 (y R a) 2 (1 - r 2 ) 3 

y R a — > 0 



(19) 



Finally, based on the rela- 
tionships (18) and (19), we 
conclude that 



Choose a Cartesian coordi- 
nate system OxiJjz, and five 
point-source locations, name- 
ly 0(0, 0, 0) , A^l, 0, 0) , A,(0, 

1, 0) , and A 4 (Q, 0, 21 ) , which 
are at (unknown) distances 
r 0 , r ± , r 2 , r 3 and r 4 , respe- 
ctively from the sphere's 
center K. The parameter 1 is 
a chosen fixed length. The 
sizes of the resulting five 

measurements E^ r .(r. |p„) are 



M . = e sc fr I 6 ) 

■» (1-mAk) 3' Pa) 



( 21 ) 



Therefore the measure of the an d 
scattered field E 3C (r 0 |p t ) is 



r, rr 

o = — b = — 

+ 7 l' \i 



( 22 ) 



Where j = 0, 1, 2, 3, 4 . Therefo- 
re we have the following five 
measurements 



b Pj 

M = — — 

J ( p\-b 4 ) 3 



(23) 



with Pj>b 2 >Q. Applying the 

law of cosines in triangle 
K0A 4 and using that 



n = pA 



(24) 
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we obtain 

p 4 =2 + 2p 3 -p a (25) 

Because of (24), we take 

p] - 3pjb 4 + 3 p 2 b s - 

,12 ,b Pi ( 26 ) 

- h 12 - b 5 ^ = 0 

Thus in equation (26) iden- 
tified the p, with j = 0, 1, 2, 

3, 4 . If the radius a of the 
sphere is known, solve the 
system of six algebraic equa- 
tions (25) and (26). Also if 
you consider that the radius 
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a is too small so that a«l, 

then the relationship (23) 
shows [3], 

h 5 

M. = -=+- (27) 



so we have the following sy- 
stem of equations 



p 4 = 2 + 2p 3 - p 0 
P 5 jMj - b 5 = 0 



where j = 0, 1, 2, 3, 4 . 
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